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Abstract 

The problem of describing the behavior of the solutions to the Navier- 
Stokes equations in three space dimensions has always been borderline. 
From one side, due to the viscosity term, smooth data seem to produce 
solutions with an everlasting regular behavior. On the other hand, the 
lack of a convincing theoretical analysis suggests the existence of possi- 
ble counterexamples. In particular, one cannot exclude the blowing up 
of solutions in finite time even in presence of smooth data. Here we give 
examples of explicit solutions of the non-homogeneous equations. These 
are defined on a Hill's type vortex where the flow is rotating and swirling 
at the same time, inducing the flux to spiraling at a central node. Despite 
the appearance, the solution still remains very regular at the agglomera- 
tion point. The analysis may lead to a better understanding of the subtle 
problem of characterizing the solution space of the 3D Navier-Stokes equa- 
tions. For instance, this result makes more narrow the path to the search 
of counterexamples built on the stretching and twisting of vortex tubes. 
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1 An instructive set of explicit solutions 

The literature on fluid dynamics offers a wide range of exact solutions, suitable 
to almost any kind of applications (see for instance |13) for a review). Many 
examples are built starting from vector spherical harmonic functions, which 
are generally used as expansion basis for more advanced scenarios (see, e.g., 
[T^). Here we deal with some of these functions, which have not to be confused 
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with their simpler version related to the scalar Laplace equation in spherical 
coordinates (see for instance [2], p. 121). For this reason, this section is going 
to be a bit technical. We discuss a specific situation from the fluid dynamics 
viewpoint, in order to prepare the ground to generalizations that will allow in 
section 3 to construct families of solutions defined on evolving vortex rings. 

First of all, we introduce the family of Bessel functions of the first kind: 

m + ^ - - -Mr) r>0 (1.1) 

where, in our context, a is a real positive parameter. It is known that, for r 
tending to zero, the following asymptotic estimate holds: 

^W-wT^Q" (1-2) 



r(a + 1) V2 

where F is the Gamma function. Moreover, Jq, turns out to be bounded, since it 
decays to zero as ^J'^pKr^ for r — > +c», with an oscillating behavior. For other 
properties on Bessel functions the reader is addressed to [14] . 

Successively, it is necessary to introduce the set of Legendre polynomials: 
{\-x^)Pl[{x) - '2.xPl^{x) + n(n + l)P„(a;) = a; G [-1, 1] n>0 (1.3) 
By the substitution x = cos the above equation takes the form: 

-l-^fsin^^ P;(cos0)) = n(n + l)P„(cos0) (1.4) 
sm d ou \ / 

We are now ready to work in spherical coordinates: 

{x, y, z) — (r sin cos (/i, rsinflsin^, rcosO) (1-5) 

with < (/)< 27r, < 6* < TT and r > 0. 

Let us start by fixing n > 1 and defining the stationary vector field w = 
{wi,W2,W3) whose components are: 

w\ = n(n + l) -= — _P„(cos6') 

W2 = (4+1/2 W + Yr J Pni^O^Q) 

^^3 = (1.6) 

Note that there is no dependence on the variable 0. Moreover, there is cylindrical 
symmetry along the z-axis. Later, in section 4, we will generalize this setting. 

Afterwards, one defines the time-dependent field v = (wi,t'2,'53) given by: 
V = e"''*w t > (1.7) 
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The constant i/ > will play the role of viscosity parameter. 

We start by discussing the regularity of w. Of course, the same arguments 
wiU apply to v, Vt > 0. It should be clear that w e C°°(R3 - {0}), Vn > 1. 
Less evident is the situation at the origin (r — 0). For n > 2, thanks to (|1.2p we 
deduce that w tends to zero. Therefore, a continuous prolongation is obtained 
by setting w = at the origin. The differentiability follows from an analysis 
similar to that of the more delicate case n — 1, to which we are now going to 
concentrate our attention. 

Let us first note that Pi (cos 0) — cos 6 and P{{cosd) — 1. Moreover, it is 
known that: 

[2 sinr - rcosr 
J3/2{r) = W j= (1.8) 

Thus, in a neighborhood of the origin one gets: 



ir^Mo + ^l -?>/lfi- 4+1^1 (1.10) 



The second component of w follows a similar asymptotic behavior 

/f \^'^/'^^ ' ■ 2r J 3 V TT V' 5 ■ 280 

The important fact is that the coefficients in (|1.9p and (|1.10l) not depending on 
r are the same. 

In order to understand what is actually happening for r = it is better to 
argue in Cartesian coordinates. The local orthogonal triplet of the spherical 
reference framework, has the expression: 

r = (sin cos 0, sin sin 0, cos 9) 

9 = (cos 9 cos 4>, cos 9 sin 0, — sin 9) 

4> = (— sin0, COS0, 0) (1-11) 

with 



sin 9 = \/ x'^ + j \/ + + cos9 = zj\J + + 

sin fp = y / x'^ + y'^ cos cj) = xj \/ x"^ + 7/^ (1-12) 



Neglecting higher order infinitesimal terms, near the origin one has from (|1.9p 
and ([LTU)) : 

w = ^ir + zi^e + zia'/' « ^, 1-^-^^") (1.13) 

^ ^ 3 V TT VlO 10' 10 5 j ^ ' 

By placing at the origin the vector (0, 0, 2\/2j'i^/W) (written in Cartesian coor- 
dinates) one gets a continuous vector field. A more careful analysis shows that 
w can be differentiated infinite times. A plot of w near the origin of the plane 
(/) = is shown in figure 1. Due to the decay at infinity of Bessel functions, the 
whole field is bounded in R"^. 
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Figure 1: Qualitative displacement of the field w near the origin for n = 1. The plot 
shows the projection on the plane (a;, 0, z) obtained from (|1.5p by setting = 0. 



The vector fields previously introduced satisfy a series of striking results that 
we are going to recall. The first important relation is that, independently of n, 
V has zero divergence: 



divv = 



V t > 



The proof is obtained by a direct check. Indeed, for a — n 
(11.41) one recovers > and r > 0: 



(1.14) 
1/2, thanks to 



divv 

n(n + 1)9 



j,2 



1 



d 



r*2 Qj* 

1 



r sin B dO 
iV^Jc) P„(cos(?) 



(sin 9 V2) 



J' 



J OL 

27 



d_ 



(sin^e P^(cos( 



(1.15) 



Equation (I1.14p is also true at the origin. To check this, one can use either (|1.13p 
or symmetry arguments (from the physical viewpoint in any neighborhood of 
the origin the inbound flow equates that directed outbound). 

Another important fact is that v satisfies the parabolic equation: 



i/Av = 

dt 



V i > 



(1.16) 
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The statement follows from writing the Laplace operator in spherical coordi- 
nates. Considering that there is no dependence on the variable </>, the expression 
of the operator looks rather simplified: 



(Av) 



1 



d_ 

dr 



,dvi 
dr 



1 d ( . dv^ 
sine 86 V'"" 89 



2 «i + ^ 



8V2 



COS ( 



V2- 



(Av) 



d_ 

dr 



,dv2 
dr 

2— 

de 



\^d9 



sm ( 



^ dv2 
09 



V2 

sm^9 



1 dv2 
V dt 



(Av)3 



1 dv3 
V dt 







(1.18) 



(1.19) 



The check of the above formulas, by direct substitution of the expressions of vi 
and V2, requires however a generous amount of patience. We do not report such 
details here. We just recall that the equations (|1.1|) and (|1.3|) must be used with 
profusion. A simpler verification is obtained by recovering a potential vector A 
such that w = rotA (see section 3) and showing that A = rotw. This implies 
that —Aw — rot(rotw) = w (recall that w has zero divergence), and this last 
relation implies (|1.16p . 

Similar solutions can be found in the 2D framework. For instance, working 
in polar coordinates (r, 9), we define for n > 1: 



(wi,W2) = (n^^^^TJ COS 9), - J'(r) sin9 T'(cos9) 
\ r 



(1.20) 



Again, one has divw = and 



where Tn is the n-th Chebyshev polynomial. 
—Aw = w. 

After having introduced the velocity field v, we need to explain its role in the 
framework of the Navier-Stokes equations. This will be the subject of the next 
section. The origin must be regarded as the bottleneck of a generous amount of 
fluid entering for z < and exiting for z > 0. Nevertheless, this is not enough 
to generate singularities in the field or its derivatives. Such an observation 
continues to be true even by introducing a swirl in the 4> direction, as we shall 
see in section 4. 



2 Connections with the Navier-Stokes equations 

A further property of the field w is that its rotational has a very simple expres- 
sion: 

rotw = (o, 0, — =J„+i/2(r)sine' P,;(cos6')) (2.21) 
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where the third component has been computed via the formula: 



(rotw)3 = - 
r 



d dwi 



(2.22) 



Due the asymptotic expression (|1.2p . the rotational goes smoothly to zero at 
the origin as r". Considering that w is bounded, the field w x rotw also goes 
to zero at the origin. 



We now set: 



V X rotv = — e w X rotw 



j>=- i|v|2 = - 16-2"^* |w|2 (2.23) 
By applying well-known notions of differential calculus, one may write: 

(v-V)v = jVjvp - V X rotv = - Vp + f (2.24) 
Therefore, using (|1.16p . we finally have: 
(9v 

zyAv + (v-V)v ==-Vp + f Vt>0 (2.25) 

at 

Recalling the divergence-free condition (|1.2p . we may argue that the couple v, 
■p satisfies the 3D Navier-Stokes problem with the forcing term f and the initial 
condition v = w at time t = 0. Note that our pressure is negative, but this is 
not a problem, since a scalar potential is defined up to an additive constant and 
|v| is bounded in the whole R'^ space. 

The choice of p looks physically correct. The quantity ipjvp, where p is the 
mass density, is related to kinetic energy. In the incompressible case, which is 
the one we are studying here, p is constant in R'^ (assume p = 1). The setting is 
in agreement with the Bernoulli principle, stating that the quantity ipjvp -V v 
is preserved along the stream-lines. Thus, the balance between — Vp and the 
gradient part of (v • V)v only involves energy properties belonging to the fiuid, 
while the extra-force f (not of conservative type in our case) represents external 
constraints. 

Including potential functions of conservative fields, such as gravity, in the 
pressure term is an allowed procedure (in this fashion, when considering the case 
of gravitational potentials, pressure directly turns out to decay with altitude). 
Effective forcing terms can be built with fields having the rotational different 
from zero. Note for instance that, if ()2.23|) was replaced by the setting f = 
and (v- V)v — — Vp, we would get an incompatible situation since the advective 
term is not the gradient of any potential. 

Our solution is defined in R'^. Of course, a fictitious boundary- value problem 
on an open domain il may be obtained by assigning on the restriction of the 
field V, Vt > 0. When is a sphere of radius R centered at the origin, there are 
two interesting situations to be mentioned. If i? is a zero of Jq (a = ri + 1/2), 



6 



we have that the radial component on dfl is identicahy zero, so the boundary 
field is tangential. The boundary datum is purely radial if instead J? is a zero 
of + Ja/2r. 

Assuming that R is the first zero of Jn+1/2 (the most interesting case is 
n = 1, where R is approximately equal to 4.49), we can interpret the solution as 
a stable vortex ring fl. Tangential velocity is imposed on the boundary and the 
term f plays the role of centripetal force keeping the fluid in stationary circular 
motion. In truth the domain is a degenerate toroid where the central hole is 
reduced to a segment of the vertical axis z, i.e., we are simulating a Hill's type 
vortex (see, e.g., [I], p. 175). 

The same construction can be considered in the 2D case (see (|1.20p ). The 
corresponding f is given by J„(r) sin6' Tl^{cos6){—W2,wi), which is tending to 
zero at the origin (see (|1.2p ) for any n > 1. 

Let us also note that in order to get a time decay at infinity for the solution 
proportional to e"''*, as the dissipative term prescribes, our forcing field f must 
go to zero at least as e"^''*, i.e. quadratically faster. This observation indicates 
that the time smoothness properties of the right-hand side should enter more 
formally in the a priori estimates, if one looks for regularity theorems. 

In the theoretical analysis of the Navier-Stokes equations, the gradient terms 
are usually eliminated by projection into the space of functions with null diver- 
gence. An alternative approach is to differentiate the momentum equation. By 
taking the rotational on both sides, one gets for some right-hand side f: 

— i^Auj + (v • V)u; - (lj ■ V)v = rotf (2.26) 

dt 

where the vorticity a; automatically satisfies diva; = 0. In this way we again 
eliminate the gradient term V(^|vp +p)- Note that u) does not satisfy homo- 
geneous boundary conditions, making difficult the theoretical analysis of (|2.26p . 

From the explicit expression of w, a certain number of collateral problems 
can be examined. For instance, one can try to find v and p such that: 

d\r 

— - i^Av + (v • V)v = - + f > 

divv = V<> and V = vo = for < = (2.27) 

Here f = — 7W x rotw, where w is defined in (|I.6|) and 7 is a given positive 
function of time with 7(0) = 0. We suppose that 7 grows for a while, then 
it decays fast to zero. Both f and Vq are defined in R"^. The solution v is 
expected to develop the same vortices as the field v introduced before. This 
should happen as soon as i > 0. Things may however depend on the rate of 
growth of 7. We guess (but this is just a conjecture to be verified with numerical 
simulations) that, if 7 follows a slow growth at the beginning, viscosity might 
dominate smoothing out the solution. Alternatively, if 7 has a rapid increase, 
the rotatory behavior may catch up and persist. 
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A variant of the example above is reformulated by replacing f in (|2.27p by 
a uniformly approximating sequence fr of functions. We can let the parameter 
T depend on t, in such a way that the limit f ~ — w x rotw is reached in a 
finite time T. It is possible to define f^- in the time interval ]0,T[ in order to 
get solutions representing classical vortex rings (i.e., closed vortex tubes shaped 
like a doughnut) tending to a Hill's type vortex as i — )• T. We study in the next 
section a general recipe to generate this kind of sequences. 

3 Solutions on evolving rings 

Let us now introduce a general approach to get families of solutions defined on 
rings. We first observe that a vector potential A may be assigned to the field 
w in (|1.6p . This is such that w ~ rotA and explicitly given by A = (0, 0, A3), 
where: 

A3 = — J„+i/2(r)sin0 P;,(cos^) (3.28) 

Information on vector potentials can be for instance retrieved in [T, p. 53, where 
connections with magnetic fields are also mentioned. It turns out that A is zero 
along the z-axis. 

Comparing with (|2.2ip . we deduce that A = rotw. As a consequence, by 
observing that rot(rotA) = — AA + V(divA) = — AA, the following relation 
holds: 

- AA = A (3.29) 

Therefore A is an eigenfunction of the vector Laplacian, corresponding to the 
eigenvalue A = 1. More precisely, the sphere centered at the origin of radius 
R (the first zero of J„+i/2), deprived of the vertical axis, is such that the small- 
est eigenvalue of the operator —A is equal to 1. The homogeneous boundary 
conditions associated to such an operator are of Dirichlet type, including in the 
boundary also the vertical segment at the intersection of the sphere and the 
z-axis. They can be however substituted by tangential type conditions: 

rotA • n = (3.30) 

with n denoting the outward normal field on dfl, that is what we got for the 
field w. 

In order to examine other similar situations, we work with cylindrical co- 
ordinates {r,z,(j)). We look for vector potentials of the form A = (0,0,^3), 
where A3 does not depend on (p. In addition we would like to impose p.29p 
and p.30|) . We are then requiring that A = 1 must be the first eigenvalue, so 
restricting a bit the choice of the domain fl. In terms of the only unknown A3, 
the vector differential problem becomes: 

( dr"^ ~^ r dr ^ dz^ ) ^ «^ ) 
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One can easily check that the two other components of AA are zero. The scalar 
equation p.31|) is either subject to homogeneous Dirichlet boundary conditions 
or to the boundary constraint: 



The above equations can be simplified further by setting B — rA^. By this 
substitution our problem takes the form; 



where S is the section of the toroid region il. Note that B is determined up to 
multiplicative constants (related to the intensity of the final velocity field). On 
the left-hand side of p. 331) we have a positive definite operator, as a result of 
little manipulation after multiplying by B/r and integrating on S. 

Since there is no dependence on (j), we can simulate flows on vortex rings 
having a symmetry axis coincident with the z-axis. An analysis of this kind 
has been considered in [5] from a different perspective, where inviscid fluids are 
coupled with electromagnetism. Note that the example introduced in section 1 
for n = 1 is related to famous Hertz solution, aimed to describe the electromag- 
netic field emanated by an infinitesimal dipole. Extensions, examined through a 
series of numerical tests, have been successively investigated in [3]. The results 
and the techniques there presented may be suitably adapted to the study of 
vortex formation and stability. 

For example, let S be the circle of radius R, situated in the plane (r, z) and 
centered at a point of the semi-axis r > very far away from the origin. By 
neglecting the terms in p.31|) containing r at the denominator, the function A3 
solves the 2D Laplace eigenvalue equation. The minimum eigenvalue in S is 
approximately equal to (^/i?)^ (where ^ « 2.4 is the first zero of Jq), so that, by 
taking R « 2.4, we get a reasonably good solution of problem p.29p . Indeed, 
the 3D ring with circular section S and vertical axis z is a domain where the 
equation — AA = A can be almost perfectly solved in terms of orthogonal basis 
functions. There are infinite other rings, more or less of the same shape and size 
as il, where (|3.29l) holds once one knows the eigenf unctions of (I3.33|) . When the 
section S approaches the origin, we can fiatten the part of it facing the z-axis, 
while preserving the minimum eigenvalue A = 1 (see figure 2). The result is the 
transformation of standard ring into a Hill's type one. 

The flow in ft may be defined by setting v = e^'^'rotA. Thus, as done 
in section 2, the full set of Navier-Stokes equations (including the divv = 
condition) is satisfied by writing: 

— - i^Av + (v • V)v = |V|vp - V X rotv = - Vp + f (3.34) 

where the gradient term on the right-hand side has been assimilated to a gradient 
of pressure, while the rest is a forcing source. This last centripetal acceleration 




(3.32) 





(3.33) 
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is the one needed to maintain the flow constrained in the ring fi. The energy of 
the flow is exponentially fading in time due to viscosity. In the above mentioned 
procedure it is not clear how to prolong the field v outside i}. Realistically, the 
external fluid is dragged by viscous effects. The fluid across the boundary dfl 
is continuous but may lose regularity there. This is the same difficulty one may 
find in the study of typhoons, at the interface between the central core and the 
peripheral dragged flow. 

The request that the eigenvalue A remains equal to a given constant is quite 
important in the framework of electromagnetic applications (where A = and 
c is the speed of light), but it is not crucial in fluid dynamics. We just have 
to pay attention when passing to the parabolic equation (dv/dt) = i/Av. For 
any flxed t, let w be the stationary solution satisfying: —Aw = A(t)w. As 
previously done, one sets v = 7(i)w. The fimction 7 has to be then determined 
according to the relation: 

^ = -I/A7 for te]0,T[ (3.35) 

yielding the exponential 7(t) = e""^* in the particular case A = 1. 



Figure 2: Sections of evolving tori. Starting from a classical ring having circular sec- 
tion (right), one can deform the shape and obtain through a sequence of intermediate 
rings a spherical doughnut (left) where the hole is reduced to a vertical segment. 



A classical way of looking for troubled solutions of the Navier-Stokes equa- 
tions is to let evolve regular vortex tubes in order to create a singularity at some 
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point. Considerations about this hypothesis are reported for instance in [B], [S], 
[9], [To] and [11], section 5.4. Some theoretical imphcations are examined in [4]. 

Therefore, from our viewpoint, one can imagine a smooth ring driven by an 
internal force f^- changing with time. In a finite time, the process ends up with 
the flow constricted into a Hill's vortex, where the set of equations are actually 
solved on a spherical domain. Note that, due to incompressibility, the volume 
of these rings must remain constant. This is also true because we are imposing 
tangential boundary conditions (see (I3.30|) '). therefore the fluid cannot escape. 

In the plane (r, z), a possible transition from a classical vortex ring displaying 
circular section into a spherical vortex is shown in figure 2. The evolution of 
the boundary may be described for t G [0, 1] by the formula: 

[V{r-s{t))^ + z^y'' = {r-s{t))^-' (3.36) 

where s > is a time-dependent decreasing function satisfying s(l)=:0. For 
t = 0, we have the circumference of radius ^ centered at (i + s(0),0). After 
a sudden change of topology, for t = 1 wc obtain the semi-circumference of 
radius 1 centered at the origin. This construction is only given to qualitatively 
illustrate the situation. We do not know for example if the family of rings built 
by rotating the domains of figure 2 around the z-axis preserves the eigenvalue 
A = 1 as prescribed by p.29p . 

It seems however that no irregularity turns out to be present at the end 
of this process, since the solution presented in section 2 is smooth up to the 
origin. This does not prove that any attempt to collapse a ring hole into a point 
(or a segment), according to the equations, is going to fail. Nevertheless, the 
percentage of success in the search of anomalous solutions of this type becomes 
very low. 



4 A Hill's type vortex ring with twist 

We finish by retouching the solution defined in section 2. We now add a com- 
ponent along the (f) direction. Thus, we go back to spherical coordinates (r, 0, 0) 
and recall the field w defined in (|1.6p . For simplicity we only discuss the case 
n = 1. We now define: 

J3/2{r) . \ 
wi, W2, cr — — — smO \ (4.37) 



where cr is a given number. Practically, the stationary flow, while behaving as 
in figure 1, also moves around the vertical axis. The third component in (|4.37p 
is zero for r = by virtue of (|1.2p . By computing the rotational, we now get: 

A = curlw= ('^^:^cos0,-^U/,(r) + :%^)sin0, (4.38) 
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where A3 is given in p.28p . It turns out that w = curlA, so we are deahng 
with the same kind of solutions previously studied. Such solutions are defined 
everywhere. However, by denoting with R the first zero of J3/2 (which is ap- 
proximately equal to 4.49), the solution can be restricted to the sphere centered 
at the origin of radius R. 

Though w still does not depend on (f), the modified fiow field twists in the 
direction of the vector cf) (see figure 3). Note that the boundary conditions to 
be assigned on the surface of the sphere of radius R only contain the component 
along 9. As figure 3 indicates, the field has only a vertical (along the z-axis) 
component at the origin. In the immediate surroundings, a component along 
(p develops, which can be very strong in magnitude since its strength depends 
on cr, which is arbitrary. Therefore, sharp layers are observed near the origin, 
although the solution continues to be smooth. Together with the origin, other 
critical regions of sharp layer formation are at the poles of the sphere. 

The trajectory of a particle is practically untwisted near the external bound- 
ary. In this situation the orbits are mildly precessing. Nevertheless, as the 
particle gets to the pole, is pushed towards the center and starts spiraling until 
it reemerges near the opposite pole. For particles placed inside the sphere far 
from the axis and the boundary, the orbits tend to follow a helicoid with cir- 
cular section, turning around the z-axis. There might be no periodicity in this 
complex motion. 




Figure 3: The field of the new spherical vortex combines the component describing 
the rotation within each section with a perpendicular one following closed circular 
patterns around the vertical axis. The horizontal velocities are different depending 
on the distance from the origin. In the neighborhood of the origin (picture on the 
right), the trajectories pass with continuity from a pure vertical motion to spirals. 
The variation can be particularly significant, especially when a is large. 
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When we compute the product w x rotw, it should be evident that its first 
two components tend to zero for r — )• 0. The third one takes the form: 

(w X rotw)3 = W1A2 — W2A1 (4.39) 

A careful analysis shows that this also tends to zero for r — )■ 0. We can argue 
locally in the neighborhood of the origin by working in Cartesian coordinates. 
In this way, excluding third-order terms, the expression given in (|1.13|) has to 
be modified as follows: 

Considering that rotw « {2\/2/ 'i^){~y /2^x /2,2a), one ends up with the ap- 
proximation: 

8 / (4(7^ — l)x ayz xz"^ ^ , 2 2\ 



-i/Aw - w X rotw w — -^^ — + -2- H 1 (z^ -I- w^), 

97r V 2 5 20 10^ ' 

(4(7^ — l)y axz yz"^ V / 2 2\ /- 2 2\ i 



5 ■ 20 ' 10(^ 20^^ +^^) + ^2V2j (^-^^^ 

The special choice a = 1/2, that allows the above terms to tend to zero faster at 
the origin, worth to be recorded. For this value of tr, the divergence of w x rotw 
is zero at the origin. 

An interesting collateral problem to be studied is when the initial solution 
of the homogeneous (f — 0) Navier-Stokes equations it taken to be equal to w. 
In addition, the fiow is constantly supplied at the boundary of the sphere of 
radius R with the stationary values of w. In absence of the forcing term it is 
reasonable to expect a decay to zero at the origin due to viscosity. This should 
be true independently of a, but strange surprises could also happen. It would 
be wise to run some numerical tests to verify the behavior. 

Of course, one can add a twisting also to the annular flows introduced in 
section 3. In cylindrical coordinates it is enough to search for solutions of 
the form w = (0,0,11)3), where ws does not depend on (p. If we require that 
simultaneously A = rotw and w = rot A, we arrive at the equation: 



d^W3 1 dw3 W3 d^ws 

Qj.2 J. Qj, J.2 Q^2 



W3 (4.42) 



to be solved in the section S. The solution is determined up to a constant 
a. By imposing homogeneous Dirichlet boundary conditions we get that the 
flow is swirling only at the interior of the ring. For a ring with small circular 
section with respect to the major diameter, the solution of the full Navier- 
Stokes equations leads to a component along the </> direction similar to that of 
a Poiseuille type flow. As before, the section S may be deformed depending on 
time as suggested by the plots of flgure 2, so ending up with a configuration as 
the one depicted in figure 3. 
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In conclusion, we proposed a set of explicit solutions of the non-homogeneous 
Navier-Stokes equations expressed in spherical coordinates with the help of 
Bessel functions. These are defined in the whole 3D space, but they can be 
easily restated in a spherical domain subject to suitable tangential conditions. 
The flow rotates forming a Hill's vortex remaining smooth up to the inner point. 
A twist may be also introduced, but this does not alter the regularity of the 
solution. This example suggests that the search of blowing up solutions of the 
3D Navier-Stokes equations is not going to be easy, since most of the counterex- 
amples are built on the idea of crashing vortex tubes in order to detect possible 
degenerations at some meeting point. Together with the mentioned examples, 
we also proposed a way to explicitly construct vortex rings, that can be helpful 
in various other circumstances, both for theoretical and computational aspects. 
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